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SUMMARY 


The Federal Aviation Administration (FAA) has Implemented an automated, time- 
based metering form of air traffic control with profile descent procedures for 
arrivals into the terminal area. These concepts provide fuel savings by 
matching the airplane-arrival flow to the airport acceptance rate through 
time— control computations and by allowing the pilot to descent at his 
discretion from cruise altitude to a designated metering fix in an Idle- 
thrust, clean-configuration (landing gear up, flaps zero, and speed brakes 
retracted). Substantial fuel savings have resulted from these procedures, but 
traffic control (ATC) workload is high since the radar controller 
maintains time management for each airplane through either speed control or 
path stretching with radar vectors. Pilot workload is also high since the 
pilot must plan for an idle— thrust descent to the metering fix using various 
rules of thumb. 

The National Aeronautics and Space Administration (NASA) has developed an 
airborne descent algorithm compatible with time— based metering procedures and 
profile descent procedures designed to improve the accuracy of delivering an 
airplane to a metering fix at a time designated by the ATC system. This 
algorithm provides open-loop guidance for an airplane to make an idle-thrust, 
clean-configured descent to arrive at the metering fix at a predetermined 
time, altitude, and airspeed. The algorithm may also be used -for planning 
fuel-conservative descents when time is not a consideration. 

The algorithm was programmed on a small programmable calculator for use with a 
DC— 10 airplane. This report describes the descent algorithm computations and 
the required vertical performance modeling for the DC-10 airplane. 







INTRODUCTION 


In an effort to Improve the efficiency of terminal area operations, the 
Federal Aviation Administration (FAA) has Implemented an automated tlme~based 
metering form of air traffic control with profile descent procedures* The 
time-based metering concept is based upon airplane arrivals crossing a 
metering fix (typically 30 to 40 n» mi. from the airport) at a specified 
altitude, airspeed and time. This results in derandomizing the arrivals at 
high altitudes resulting in a reduction of the low altitude, high fuel con- 
sumption flight required for sequencing airplanes to a common final approach 
P®bh. In addition, time- based metering allows the traffic to absorb delays at 
cruise altitudes or while on the ground prior to takeoff, resulting in even 
greater fuel savings. The proper sequencing and spacing of enroute traffic 
also allows for increased airport productivity (ref. 1 and 2). The profile 
descent procedure allows the pilot to plan and fly a descent which is fuel 
conservative for his particular airplane resulting in an additional fuel 
savings . 

With the current time-based metering/profile descent procedures, the air 
traffic controller is responsible for the time management of each aircraft. 

The controller can adjust the airplane's time of arrival at the metering fix 
by increasing the flight path through heading changes or by requesting the 
pilot to change speed. The pilot is responsible for crossing the metering fix 
at the proper airspeed and altitude and must plan the descent carefully if 
fuel is to be conserved. The time management and the descent planning are 
both workload intensive and interdependent. Since limited or no guidance is 
available to either the controller or the pilot, their tasks must be 
accomplished independently through various rules of thumb and past experience. 
With this present operational concept, airplanes will typically cross the 
metering fix with a time accuracy between 1 and 2 minutes (ref. 3). 

the summer of 1979, the National Aeronautics and Space Administration 
(NASA) developed and flight tested in its Transport Systems Research Vehicle 
(TSRV), previously designated the Terminal Configured Vehicle B-737 research 
airplane, a flight -management descent algorithm designed to provide closed- 
loop guidance for a fuel conservative descent and to reduce the metering fix 
crossing time dispersion. The flight management algorithm descent compu- 
tations were based on an idle-thrust, clean-configured descent (landing gear 
up, flaps and spoilers retracted), to arrive at the metering fix at the proper 
altitude, airspeed, and ATC— designated time. The results of the flight tests 
showed that the closed-loop guidance provided by the guidance and display 
system could reduce the crossing time dispersion to approximately 12 seconds 
(ref. 4). 

This research was continued in June of 1981 with a T-39A (Sabreliner) airplane 
to determine if similar results could be obtained with open-loop guidance and 
a conventional complement of cockpit instrumentation (ref. 5). A version of 
the flight management descent algorithm was implemented on a small program- 
mable calculator. Open-loop guidance was provided by the calculator in the 
form of the Mach number and airspeed at which the descent should be flown and 
the point at which the pilot was to reduce the thrust to flight idle and begin 
the descent. The descent was then flown with reference to the airplane Mach 


3 



and airspeed Indicators. Flight tests using this open-loop guidance resulted 
in a time dispersion crossing the metering fix of approximately 20 seconds. 

Having determined the viability of the open-loop descent guidance, additional 
research has been planned to evaluate the feasibility of such a descent 
planning tool in an airline operational environment. The airplane selected 
for this study was the McDonnell Douglas DC-10. This report contains a 
description of the programmable calculator software and of the DC-10 descent 
performance model used in the algorithm computations. A later report will 
present the results of the feasibility study. 

Use of company names or designations in this report does not constitute an 
official endorsement of such companies or products, either expressed or 
implied, by the National Aeronautics and Space Administration. 
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SYMBOLS AND ABBREVIATIONS 


ATC 


^0 

^1 


air traffic control 

coefficients for quadratic curve fits of altitude as a 
function of time 

modeled vertical speed at sea level, ft/sec 
slope of the linear model for h^^g, ft/sec/ft 


constant' in the model for hj^ , sec^/ft 
constant in the model for h^ , ft 


DME 


^w,h 


w,s 


distance measuring equipment 

magnetic wind direction evaluated at altitude h, deg 

magnetic wind direction computed for sea-level altitude, 
deg 


dD^ wind direction gradient with respect to altitude, deg/ft 

dh 


dh 

GW 


GSc 

H 

h 

^c 

^XO 


ft 


wind speed gradient with respect to altitude, knots /ft 


gross weight , lb 

ground speed at cruise altitude, knots 
pressure altitude, ft 
geopotential altitude, ft 
cruise altitude, ft 
metering-fix altitude, ft 

altitude at transition from constant -Mach descent to 
constant-airspeed descent, ft 

rate of change of altitude, ft/sec 
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^lAS, 

d 

IAS 

lAS^ 

Initial 


IAS 


d,i 


IAS 


IAS 


*d,max 

d,min 


^ASmf 
IDL DME 


K 


gw 


M 

M/IAS 

initial 

MSL 

“d 

^d, Initial 
OAT 


w,h 


rate of change of altitude evaluated at indicated airspeed 
lAS^, ft/sec 

rate of change of altitude evaluated at Mach number M, , 
ft/sec 

indicated airspeed, knots 

indicated airspeed used during descent, knots 

initial descent indicated airspeed for speed iteration 
computations , knots 

descent indicated airspeed computed on ith iteration, knots 

maximum operational descent indicated airspeed, knots 

minimum operational descent indicated airspeed, knots 

indicated airspeed to cross metering fix, knots 

calculator display showing point where thrust should be 
reduced to flight idle, n. mi. 

interpolation factor computed for speed iteration purposes 
gross-weight multiplication factor for altitude rate 
distance between entry fix and metering fix, n. mi. 

Mach number 

Mach number and indicated airspeed 

initial Mach number and indicated airspeed for speed 
iteration computations, knots 

mean sea level 

cruise Mach number 

descent Mach number 

initial descent Mach number for speed iteration purposes 
outside air temperature, 

wind speed evaluated at altitude h, knots 
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w,s 

t 

TAS 

TRK 

TSRV 

Tc 

"*^i8a,h 


T» 

o 


st,c 

”^st,h 

*^E 

*^EF 

VOR 


W, 


H,c 


W, 


H,h 


AA, 


AT 


At 

At 


j 

req 


wind speed computed for seal“level altitude, knots 
time, sec 

true airspeed, knots 

airplane magnetic track angle along ground, deg 
Transport Systems Research Vehicle 
temperature measured at cruise altitude, 

International Standard Atmospheric temperature at altitude 
h, °K 

standard sea-level air temperature, °K 

nonstandard sea-level air temperature, °K 

static air temperature measured at cruise altitude, °K 

static air temperature at altitude h, °K 

time error for descent -speed convergence criteria, sec 

time that entry fix was crossed, hr;min:sec 

very high frequency omnirange navigation radio 

difference between actual and computed ground speeds at 
cruise altitude, knots 

head— wind component along airplane ground track evaluated 
at altitude h, knots 

acceleration, knots /sec 

length of path segment j , n. mi . 

difference between actual temperature and standard 
temperature , °K 

time required to fly on path segment j , sec 

time required to fly between entry fix and metering fix, 
sec 
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DESCRIPTION OF FLIGHT-MANAGEMENT DESCENT ALGORITHM 
Description of General Profile 

The flight -management descent algorithm computes the parameters required to 
describe a seven-segment cruise and descent profile (fig. 1) between an 
arbitrarily located entry fix and an ATC-defined metering fix. The descent 
profile is computed based on empirical modeling of airplane performance for an 
idle-thrust, clean-configured descent. Airplane gross weight, wind, and 
nonstandard-temperature effects are also considered in these calculations. 

Figure 1 shows the vertical-plane geometry of the path between the entry fix 
and the metering fix. Each path segment, starting at the metering fix, is 
numbered according to the order in which it is calculated by the algorithm. 

To be compatible with standard airline operating practices, the path is 
calculated based upon the descent being flown at a constant Mach number with 
transition to a constant indicated airspeed and all speed reductions made in 
level flight. 

The first segment traversed on the profile is segment 7 which begins at the 
entry fix and is flown at constant cruise altitude and Mach number. Segment 6 
is a relatively short, level-flight path segment in which the pilot reduces 
thrust to flight idle so the airplane will slow from the cruise Mach number to 
the descent Mach number. Segment 6 is eliminated if the descent and cruise 
Mach numbers are the same. Once the descent Mach number is attained, the 
constant Mach descent segment (segment 5) is started. As altitude is 
decreased along this path segment, the Indicated airspeed will Increase 
because of increasing air pressure. Segment 4 begins when the desired 
indicated airspeed is attained for descent. The descent is continued along 
this segment at the desired, constant Indicated airspeed. When the metering 
fix altitude has been reached, the airplane Is flown at a constant altitude 
along segment 3 and slowed from the descent airspeed to the designated air- 
speed over the metering fix. If the metering-fix altitude is below 10 000 
feet MSL and the descent airspeed is greater than 250 knots, segments 1 and 2 
are computed for the pilot to comply with the ATC-imposed airspeed limit of 
250 knots, or less, below 10 000 feet MSL. Segment 3 then becomes a level- 
flight segment at 10 000 feet MSL where the airspeed is reduced to 250 knots. 
The descent is then continued at 250 knots along segment 2. When the 
metering-fix altitude has been reached, the airplane is flown at a constant 
altitude along segment 1 and slowed from the descent airspeed to the desig- 
nated airspeed over the metering fix. This path segment is eliminated if the 
descent and metering-fix airspeeds are the same. 

The flight -management descent algorithm can be used in either of two modes. 

In the first mode the pilot can enter the desired M/IAS to be flown during 
descent. The descent profile is then computed based on this descent speed 
schedule without consideration of a constraint on metering-fix arrival time. 
This mode would be used when time-based metering is not being used. 

The second mode was designed for time-metered operations. In this mode. 
Instead of specifying the M/IAS descent schedule, the pilot enters the time 
that the entry fix was crossed and the metering-fix arrival time assigned by 
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ATC* The descent profile is then calculated based on a M/IAS descent 
schedule, computed through an interative process, that will closely satisfy 
the crossing time specified for the metering fix. The magnitude of the Mach 
number programmed for the descent is the same as that used for the cruise 
segment. 

During the descent profile computations in the time-metered mode, a check is 
made to insure that the descent airspeed lAS^ is within the minimum and 

maximum speed limits for the particular airplane modeled. For the DC-10 
airplane, these limits were 

220 £ lAS^ ^ 350 [knots] 


There was an additional constraint that lAS^ would not be less than the 

airspeed at which the airplane was to cross the metering fix, so that extra 
fuel would not be required as thrust is added to subsequently Increase air- 
speed. If the ATC-assigned metering-fix crossing time requires a descent 
speed less than the airplane’s minimum descent speed limit, the profile is 
computed based on the minimum allowable descent speed and a message is 
displayed to the plot to "hold** (delay) for an appropriate amount of time. A 
similar "late” message is displayed with the time error if a descent speed 
schedule greater than the maximum allowed is required. 


Logic Flow of Profile Descent Algorithm 


Figure 2 shows the general logic flow of the profile descent computations. 
Pilot inputs used to compute the profile may be entered prior to flight and 
modified, as required, prior to the descent. These parameters include cruise 
altitude and Mach number, airplane gross weight, outside air temperature, 
entry-fix and metering-fix descriptions, and the course direction to the 
metering fix. In addition to these parameters, the pilot may enter either a 
particular Mach number and indicated airspeed to be used during the descent or 
the entry-fix crossing time and the ATC-assigned metering-fix crossing time. 

If the M/IAS descent speed schedule has been entered in the calculator, the 
computations will be based on a nonmetered traffic environment. The pilot 
initiates the computations by pushing the "compute" key. The descent profile 
is then computed in a single iteration, and the display gives the point where 
thrust should be reduced to flight idle to start the descent. 


If the entry-fix crossing time and the ATC-assigned metering-fix crossing time 

have been entered in the calculator, the time required to fly between the 

fixes At will be computed and subsequent calculations will be based on a 
req 

time-metered traffic environment. Once the pilot has initiated the compu- 
tations by pushing the "compute” key, an interative process is, started to 
determine an appropriate IAS descent speed that will satisfy the time 
constraints. The Mach number used for the descent is set equal to the cruise 
Mach number in an effort to reduce computational requirements and operational 
complexity. 
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The Iterative process starts with the computation of the time to complete a 
descent Initial following descent speed schedule. 


M, . , = M 

d, Initial c 


^d. Initial 


280 knots If lAS^ < 280 
MF 

lASjjp otherwise 


The initial 280 knots Is used because It Is the approximate 

mld-polnt between the maximum and minimum allowable descent airspeeds and 
because It Is a descent speed typically used by the airlines. 


A check Is made with the following transition-time inequality to determine if 
the time-convergence criteria tg has been satisfied. (For the purposes of 
these tests, tg = 5 sec.) 



[sec] 


[sec] 


If this inequality is satisfied, the computations are complete and the idle— 
thrust descent point is displayed to the pilot. If the Inequality is not 
satisfied, the descent computations will be repeated using the operational 
airspeed limits, as follows: 
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lAS^ « 
d 


A check is then made to determine if the time criteria has been satisfied or 
if a speed greater than 350 knots or less than lAS^ would be required to 

satisfy the time constraints. If the time criteria is satisfied, the idle- 
thrust descent point is displayed to the pilot. If either the upper or lower 
airspeed limit must be violated to satisfy the time constraints, the 


350 knots if 


At < 
req 


2^ At 

j=l • 


Initial 


IAS_ otherwise 

Mr 
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appropriate speed limit will be used in the descent computations and the 
resulting time error for crossing the metering fix will be displayed to the 
pilot. 

If the time criteria has not been satisfied and neither airspeed limitation 
will be violated, a revised descent airspeed lAS^ and associated descent 

time will be computed and compared to The Iterative process will 

continue until the time convergence criteria has been satisfied. 


The computation of the revised lAS^ is graphically depicted in figure 3, 

which is a plot of the time required to fly between a specified entry fix and 
metering fix at a specified cruise Mach number over the complete lAS^ range 

of the airplane. The descent airspeed is revised through a modified linear 

interpolation of the desired At within a range of time bounded by an 

req 

initial value and a computed variable value. The initial value 

At. I was the resulting time computed with the M/IAS. 4^4 

Vinitlal d.lnitlal 

descent speed schedule on the first iteration. The variable time value 


( y^ At . j is the time computed for the 
j=i ^/i-i 




descent speed schedule 


on the last (i-1) iteration. The revised descent speed schedule M/IAS^ ^ is 
computed as follows: * 


Md =Mc 

^^d,i " ^^d, initial ^^^d,i-l " ^^d.initial^*^ 

(IAS, . , - IAS, , ) 

initial 

(l)At) -(ZAt) 

\j=l J /initial \j=l ^/i-l 


[knots ] 


where K 


i: it 

1»1 ■'/initial 


- At 


req 


(e itj) 

\J =1 ^/initial \j=l J/i -1 


and 1 is the 1th iteration. 
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The last term of the computations for lAS^ is a compensation factor for 

the difference in curvature between the plot of time required to fly between 
the entry fix and the metering fix and the straight line used in the linear 
interpolation. 


Empirical Representation of Airplane Performance Characteristics 

Computer memory limitations with the programmable calculator preclude the use 
of detailed aerodynamic and performance tables to represent the airplane for 
profile decent calculations. Instead, an empirical model of the performance 
of the DC-10 airplane was developed from flight data collected during Idle- 
thrust clean-configured descents and during level-flight speed reductions. 

A portable voice recorder, a stop watch, and conventional flight Instruments 
were used to collect data for descent-performance modeling. Altitude, speed, 
temperature, and time were recorded at altitude increments of approximately 
500 feet during the descents. Speed, temperature, and time were recorded at 
10-second intervals during constant -altitude speed changes. Gross weights 
were recorded at the beginning and end of each test run. 

Constant IAS Descent Model. - The performance model for the constant indicated 
airspeed descents consisted of linear approximations of vertical speed as a 
function of altitude for a range of airspeeds between 220 knots and 350 knots. 
The vertical speed was adjusted to compensate for variations in gross weight. 

The first step in the development of the vertical performance model for 
constant lAS^ was to approximate, for each descent, altitude as a function 

of time as a quadratic equation through a least squares curve fit analysis. 

The general form assumed for this equation was 

b = a2 t^ + aj^t + Sq [ft] 


Figure 4 shows a typical plot of the data and resulting curve fit for a 
descent flown at a constant lAS^j of 280 knots. 

Vertical speed was determined by differentiating this equation with respect to 
time. This resulted in an equation of the following form: 


h = 2 a 2 t + aj [ft/sec] 


A plot of as a function of h was then developed for each descent. 

Figure 5 shows typical data for various descent speeds. These plots indicated 
that vertical speed was approximately linear as a function of altitude and 
were modeled with an equation of the form 

li = bj h + bQ [ft/sec] 
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The slope was approximately the same for all descent speeds and equal to 

-3.5 X 10”^ ft/sec/ft. The modeled vertical speed at sea level bQ was 

corrected for gross weight variations by dividing by the gross weight 
correction factor K (to be subsequently discussed). This data varied 

gw 

exponentially as a function of airspeed. A plot of this data and the 
exponential model of bQ is shown in figure 6. The resulting model of 

vertical speed, corrected for gross weight variations of the DC-10 airplane 
was 


lilAS = -3.5 X 10“^ h - 3.07783 K [ft/sec] 

d ^ 

Constant Mach Number Descent Model, - The performance model for constant Mach 
number descents consisted of parabolic approximations of vertical speed as a 
function of altitude for a range of Mach numbers betwen 0.73 and 0*85, The 
vertical speed was adjusted to compensate for variations in gross weight* 

The procedures used to develop the descent model for constant Mach number were 
similar to those used for development of the model for constant IAS descents* 

A quadratic equation of altitude as a function of time was derived for each 
descent through a least squares curve fit analysis. These equations were 
differentiated with respect to time to determine vertical speed as a function 
of time. The vertical speed was adjusted for gross weight variations by 
dividing by the gross weight correction factor K (to be subsequently 
discussed). ^ 

The vertical speed li adjusted for gross weight variation, was then plotted 
as a function of altitude. Typical plots of this data are shown in figure 
7. The resulting curves were parabolic in nature and were modeled with an 
equation of the following form: 


[ft/sec] 

The magnitude of the coefficient Cq was subjectively selected based on the 
shape of a generic parabola that would overlay the descent data. A value of 
Cq = -1.85 was selected and resulted in the parabola shown by the dashed line 
in figure 7. 

The coefficient Cj was calculated for each descent at an altitude 

approximately 2,000 feet below cruise altitude for Cq = -1.85. The resulting 

values for Cj^ are plotted as a function of the descent Mach number in figure 
8. A linear regression analysis resulted in an equation of c^ as a function 

of the Mach number. However, the slope of this model was increased slightly 
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to Insure that Cj would be greater than the maximum cruise altitude of the 

airplane (defined by the altitude limits shown on figure 8) thus Insuring that 
an imaginary root would not be obtained from the equation for ft. The 
resulting model for Is 

Cj = 25750 + 22167 


Acceleration performance model, — Acceleration performance data were obtained 
for idle-thrust clean-configured speed reductions on level flight paths for 
typical cruise and metering fix altitudes. Indicated airspeed and time data 
were recorded during the speed reductions. The Indicated airspeeds were 
converted to true airspeeds to reduce computational requirements within the 
descent algorithm. Figure 9 shows a plot of true airspeed as a function of 
time that resulted during speed reductions at various altitudes and gross 
weights. The average slope of each of these test runs was the same. Hence, 
acceleration for the DG-10 airplane was modeled as constant and equal to 


X = -1,3 [knots /sec] 


Gross weight variation, - The effects of gross weight variations on the 
airplane's descent peformance were accounted for with a single multiplication 
factor applied to the vertical performance models developed for both constant 
airspeed and constant Mach number descents. The multiplication factor K 
is a linear, nondlmenslonal expression, ^ 

The multiplication factor was derived for the DC-10 airplane by plotting the 
vertical speeds obtained during descents (conducted at the same indicated 
airspeed) with the gross weight of the airplane. Figure 10 shows the vertical 
speed at sea level as a function of gross weight for descents conducted at 
constant airspeeds of 250 knots, 280 knots, and 340 knots. A linear curve fit 
was applied to the data points for the 280 knot descents since more descents 
were flown at this speed and since a wider range of gross weight existed in 
the data. This plot also shows that the model derived for the 280 knot 
descents may be shifted vertically (maintaining the same slope) to overlay the 
descent data obtained at the other speeds. Since the same slope could be 
maintained for all airspeeds, changes to the vertical speed due to gross 
weight variations were independent of airspeed. The plot in figure 10 was 
then nondlmensionalized by dividing the abscissa by 304,000 pounds and the 
ordinate by -30.2 feet/second (bg = -30.2 ft/sec at an airspeed of 280 knots 

and a gross weight of 304,000 lbs.) This nondimens ionalization process 
allowed the gross weight variation model derived with the 280 knot descent 
data to be expressed in a form useful for descents at any airspeed. The 
resulting multiplication factor K for gross weight variations wast 

gw 

Kg^ « -3.863133 x 10"^ GW + 2.174392369 
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Approximation of True Airspeed 


It Is necessary to determine true airspeed from both Mach number and 
calibrated airspeed, as required by the path segment, so that a headwind 
component can be added to obtain ground speed for time calculations. True 
airspeed, as a function of Mach number and static air temperature Tg^ 1 ^, was 
represented by the following equation: * 


TAS = 29.04 (Tg(.^h)^-'^M [knots] 


True airspeed, as a function of calibrated airspeed and altitude h, was 
approximated as 

h < 42,000 ft 
220 < IAS < 360 knots 


TAS = 


IAS 


-4, 


[knots ] 


1 - (0.12 X 10 )h 


Wind Modeling Technique 

A two-component linear wind model was used to represent the wind speed and the 
wind direction as functions of altitude. The coefficients of the wind model 
were computed via a linear regression analysis of winds aloft reports and 
forecasts In the descent area. Winds aloft data for the linear regression 
analysis was Inserted through the calculator keyboard In a format similar to 
that used In standard aviation winds aloft forecasts. Though wind speed and 
direction from only two altitudes were required to define a wind model, the 
pilot could choose to Insert additional wind data based on both forecasts and 
pilot reports. 


The magnitude of the wind speed and the direction of the wind defined by the 
linear wind model were computed for each segment of the profile based on the 
middle altitude of each segment. The following equations were used for these 
computations. 




h + S 


w.s 


[knots ] 


w,h 


[w) ^ ^ \s 


[deg] 


A head-wind component for each segment would be computed automatically during 
the profile computations by multiplying the wind speed with the cosine of the 
angle between the airplane ground track and the wind direction. A head-wind 
component correction factor, based on the actual winds encountered during 
cruise flight could also be added to the wind model If the pilot determined It 
was necessary. The correction factor was obtained by first computing the 
difference between the actual ground speed along the cruise segment and the 
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predicted speed based on the modeled winds and cruise Mach number. This 
component was proportional to altitude and decreased linearly to zero at sea 
level (h=0). The corrected head wind component profile 

computations was defined by the following equation: 

- TRK) + (h/h^) [knots] 

Compensation for Effects of Nonstandard Atmospheric Temperature 

Various flight instruments, including the Mach meter, the airspeed indicator, 
and the altimeter, are designed to display correct indications in a standard 
atmosphere. However, standard atmospheric conditions are rarely encountered. 
This results in slight errors in indicated altitude and speed* The profile 
descent algorithm compensates for nonstandard temperatures as they affect the 
Mach number calculations and altimeter indications. 

Nonstandard temperatures are computed by the algorithm based on a standard 
atmospheric temperature model with a bias correction based on the difference 
between the actual and the standard temperatures. The standard temperature 
model is a two-segment linear profile defined as a function of altitude. This 
temperature model uses a slope equal to a temperature lapse rate of -1.978 x 

10“^ “C/ft for altitudes below the tropopause (h ^ 36 152 ft). At higher 
altitudes (h > 36 152 ft), it was assumed that flight was being conducted 
within the tropopause where temperature remains constant with changes in 
altitude. The standard temperature model is represented mathematically as: 


Tisa ^ = 216.65 h > 36 152 ft [“K] 

T^ ^ = 216.65 + 1.978 x lO"^ (36 152 - h) h < 36 152 ft [°K] 

lsa,n — 

The following bias correction AT, representing the difference between the 
actual temperature measured at cruise altitude T^ and the standard 

temperature for cruise altitude ^ is added to the standard atmospheric 

temperature profile to completely define static temperature T^^. ^ at any 

altitude h as follows: 


AT = T - T. 

c lsa,c 

[Or] 

T ^ = T, . + AT 

st,h lsa,h 

[°K] 


The static temperature is then used for conversion of Mach number to true 
airspeed. 
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Pressure altitudes H used to define the end points of each segment are 
corrected to geopotentlal altitudes by multiplying the pressure altitude by a 
temperature ratio of nonstandard and standard sea-level temperatures (ref. 6) 
as follows : 


h = H (T'/T ) [ft] 

o o 


The standard sea-level temperature Is 288.15 ®K; the nonstandard sea- 

level temperature Is computed from the static— temperature model for h *• 0. 


Computations of Descent Path 

The point where the pilot Is to reduce power to Idle thrust to start the 
descent was defined by summing the distances required to fly segments 1 to 
6. Each segment length was determined by first computing the required time to 
traverse the segment and then multiplying by the average ground speed computed 
for the segment. Times for the level— flight segments requiring airspeed or 
Mach reductions were determined by dividing the required speed change by the 
deceleration capability of the airplane. Times for the path segments 
requiring descents were determined by integrating the reciprocal of the 
vertical speed model for that segment over the altitude change required. The 
average ground speed at which the airplane was to fly each segment was 
determined by summing the computed true airspeed and the head-wind component 
evaluated for each segment. 

The cruise segment (segment 7) at level flight and constant Mach number had no 
influence on the location of the point where idle thrust was to begin. This 
segment was significant only during the time-metered mode and was used for the 
calculations to satisfy the time constraints. Segments 1 and 2 were computed 
only if the ATC-imposed limit of 250 knots indicated airspeed for flight below 
10 000 ft MSL was applicable. The details of these calculations are presented 
in the following paragraphs. 

H/lAS transition altitude. — As the airplane descends at a constant Mach 
number, the indicated airspeed increases because of an increase in the air 
pressure. The altitude at which the desired descent airspeed is obtained is 
called the M/IAS transition altitude and defines the point at which the 
constant -Mach segment ends and the constant -IAS segment begins. The general 
equation for this transition altitude was determined by equating true airspeed 
as a function of indicated airspeed and altitude with true airspeed as a 
function of Mach number and altitude. Solving for altitude results in the 
following equation to define the altitude for transition of Mach to indicated 
airspeed 

hj^Q = 1.77675 X 10^ - 

Segment 1 .— Path segment 1 is a level-flight segment on which the airplane is 
slowed from an indicated airspeed of 250 knots to the metering-fix crossing 
speed. If the metering fix crossing speed is 250 knots, or if the ATC-imposed 


(8.90046 X 10^) + (3.42936 x 10^) 


IAS 




M 


d I 


[ft] 
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250 knot maximum airspeed limit for flight below 10 000 ft MSL is not 
applicable, this segment is not computed. 

The equations for time and length in segment 1 are 

(IAS - 250) 

At, [sec] 

X [1 - (0.12 X 10 h^] 

where x = -1.3 knots/sec, and 


Ai 


1 


■(lAS^ + 250)/2 

.1 - ( 0.12 X 10 


- W, 


H,h 


MF, 


At, 


3600 


[n. ml.] 


Segment 2. - Segment 2 is an Idle-thrust descent flown at a constant 250 knots 
from 10 000 ft MSL to the metering-fix altitude. Segment 2 Is not computed if 
the metering-fix altitude is equal to or greater than 10 000 ft MSL or if the 
descent speed lAS^ flown on segment 4 is 250 knots or less. 

The equations for time and length in segment 2 are 


At, 




MF 1 


b, + b, 


L 10 000 b, 


Vo 


+ b. 


[sec] 


where 


bj = -3.5 X 10“^ 

bf, = -3.07783 K 158681x10“^) 

1 1 cm * 


Ai2 = 


250 

(l - 0.12 X 10“^^ 1°’°00 + h^ H,(10 000 + hj^)/2 

\ ^ o 


At2 

3600 


[n. ml.] 


Segment 3. - Segment 3 is a level-flight segment on which the airplane is 
slowed from the descent speed lAS^ to the metering— fix crossing speed (or 

250 knots if segments 1 and 2 are computed). 

The equations for time and length in segment 3 are 


At 


3 


^^^MF - ^^d 
x[l - (0.12 X 10"^)h] 


where x = -1.3 knots/sec, and 


[sec] 
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[n. ml.] 


AJl, 


+ lAS^ /2 


- W„ 


Ll - (0.12 X 10"^)h 


^^3 

3600 


h^, if segments 1 and 2 are not computed 


10 000 ft, otherwise 


Segment 4. - Segment 4 is an idle-thrust descent flown at a constant Indicated 
airspeed lAS^. The descent begins at the transition altitude hj^Q and ends 

at the metering-fix altitude (or at 10 000 ft MSL, if segments 1 and 2 are 
computed ) . 


The equations for the time and length of segment 4 are 


where 


and 


At, = In 
4 


^0 

.^ 0^1 ^ 0 . 



-3.5 X 10"^ 

-3 07783 K gIASj(8. 15868 x 10 [sec] 

gw 

10 000 T*^/T^, if segments 1 and 2 are computed 
hjjp , otherwise 


IAS, 


Ai^ = 


1 - (0.12 X lO"^) 




- W, 


H,(Hjjq + h)/2 


At, 

3600 


[n. ml.] 


Segment 5 . - Segment 5 is a constant-Mach descent flown at idle-thrust power 
settings. This segment begins at cruise altitude h^ and ends when lAS^ is 

attained at the transition altitude hj^Q. 


The equations for the time and length of segment 5 are 



[sec] 


where 


Cq — —1.85 

cj » 25750 + 22167 
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and 


M 3 - [29.04(T^,_3)‘^' ^ ^ ^., 

The static temperature Tg^ 5 and the head-wind component are evaluated at 
the average altitude between the cruise and transition altitudes. 

Segment 6 .- Segment 6 is a level-flight speed change from the cruise Mach 
number to the descent Mach number. If the cruise and descent Mach numbers are 
the same, this segment is not computed. The equations for time and length of 
segment 6 are 

At, = 29.04 (T 

D St,C 


M. - M 
d c 


[sec] 


where x = -1.3 knots /sec, and 




29.04 (T ) 

St ,c 


1/2 




3600 


[n. mi.] 


Segment 7.- Segment 7 is the remaining path between the entry fix and the 
beginning of segment 6. The length of segment 7 is the difference between the 
total distance between the entry fix and metering fix and the sum of the 

distances of the remaining six segments. The length is given as follows 

6 

AA_ = Aj. “ X) ^^4 
j=l J 


Segment 7 time At^ is found by dividing the distance to be flown by the 
ground speed, as follows 


At 


7 


3600 Aly 

29.04(T^^_c>^'^ «c 



[sec] 


Input /Output Requirements 

The data required for the profile descent equations are obtained from the 
preprogrammed calculator memory and from pilot entries through the keyboard 
shown in figure 11. Though all the data necessary to compute the descent are 
entered prior to takeoff, these parameters may be updated during cruise to 
obtain more accurate results. 

The wind data is entered through the keyboard and the wind model coefficients 
are automatically computed and stored in the proper memory locations. The 
wind data, correlated to altitude, is Inserted in a data format similar to 
that found on an aviation weather forecast. To Insert the wind data, the 
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pilot must first push the key labeled followed by the key labeled **WIND.” 
The display will request the altitude for the wind speed and direction data 
with the message ”H=?FT’* The altitude should be keyed into the display and 
entered into memory by pushing the "New Entry" key. The calculator will then 
request the wind direction and speed with the message "DIR.SPD?" Wind 
direction and speed should be keyed into the display and entered into memory 
by pushing the "New Entry" key. This process will be repeated until all of 
the wind data had been Inserted in the calculator. The linear regression 
analysis will be completed after the pilot Inserts a negative altitude to 
indicate that no more wind data will be Inserted. The calculator will then 
display a "WIND IN" message. 

The wind data used by the pilot to compute the wind model could contain some 
errors since that Information Is usually based on aviation forecasts or pilot 
reports. A procedure was developed that allows the pilot to modify the wind 
model with a correction factor based on the difference between the computed 
and actual ground speeds along the cruise segment (segment 7). The computed 
ground speed used in the profile descent computations may be displayed by 
pushing the key followed by the "*GSc" key. The difference between the 
displayed ground speed and the actual ground speed represents the wind 
modeling error along the magnetic course of the airplane to the metering 
fix. If the ground speeds are different, the actual ground speed may be keyed 
into the display. Then, by pushing the "New Entry" key, the difference 
between the ground speeds is computed and stored in memory for use in subse- 
quent descent and ground speed computations. 

The operational parameters affected by ATC constraints or pilot desires, and 
not accurately known until prior to the start of descent, were designed to be 
single key inputs. To enter these data, the pilot presses the particular key 
dedicated to the parameter to be changed. After the key has been pressed, the 
display will show the name of the parameter and its current value stored in 
the calculator. Another numerical value may be keyed on the display, then 
stored in the proper memory location by simply pressing the "New Entry" key. 

If the current value shown is satisfactory, no more keyboard actions will be 
required for that parameter. 

The operational parameters may be inserted in any order, or may be changed at 
any time prior to Initiating the descent calculations. When the magnitudes of 
the parameters are satisfactory to the pilot, the profile descent computations 
are initiated by pressing the "Profile" key. Computations typically require 
less than 2 minutes for completion In the tlme-metered mode of operation and 
approximately 25 seconds In the nonmetered mode. 

The operational parameters to be entered by the flight crew through the 
keyboard as well as their symbology as presented on the keyboard and the 
display, are as follows: 
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Keyboard 

symbology 

Display 

symbology 

Operational parameter 

Me 

Me 

Cruise Mach number 

He 

He 

Cruise altitude, ft 

*GSc 

GSe 

Ground speed at cruise altitude, knots 

Md 

Md 

Descent Mach number 

lASd 

lASd 

Descent indicated airspeed, knots 

Time MF 

MF TM 

Time assigned by ATC to cross metering fix, 
hr:mln:sec 

Time EF 

EF TM 

Time to cross entry fix, hr:mln:sec 

GW 

GW 

Airplane gross weight at top of descent, lbs 

OAT 

OAT 

Static outside air temperature, “C 

Metering Fix 

H MF 

Crossing altitude at metering fix, ft 

Metering Fix 

IAS MF 

Indicated airspeed to cross metering fix, 
knots 

Metering Fix 

MF DME 

DME indication defining metering— fix 
location, n. ml. 

Entry Fix 

EF DME 

DME indication defining entry-fix location - 
this mileage must be relative to same IME 
station used to define metering fix, 
n. ml. 

Entry Fix 

MAGCRS 

Magnetic course from entry fix to metering 
fix, deg 

Entry Fix 

MAGYAR 

Magnetic variation in the descent area, deg 


If the descent speed schedule has been specified, the entry-fix and metering- 
fix crossing times must remain unassigned. If these times are specified 
through the keyboard, the proper descent speed schedule will be computed and 
stored in the correct memory location for recall by the pilot. 

When the computations are completed, the display will normally show the DME 
indication where thrust should be reduced to flight idle for the descent to 
the metering fist. If the assigned metering-fix crossing time cannot be 
attained in the time— metered mode because of airplane operational speed 
limitations, a message will be displayed indicating the amount of time 
required to delay (hold) before starting the descent or the amount of time 
that the airplane will arrive late at the metering fix. 

After the profile descent computations have been completed, the value of any 
operational parameters, including those required for input, may be displayed 
by pressing the particular designated key on the keyboard. Parameters that 
may be displayed after the descent computations, and their designated names, 
are; 
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Keyboard 

symbology 


Display 

symbology 

Operational parameter 

Md 


Md 

Descent Mach number required to satisfy 
entry-fix and meterlng-flx time 
constraints 

lASd 


lASd 

Descent Indicated airspeed required to 
satisfy entry-fix and metering-fix time 
constraints, knots 

Idle DME 


IDL IME 

DME indication showing point where thrust 
should be reduced to flight Idle, n. mi. 

Late 


LATE 

Amount of time that airplane will be late 
crossing metering fix, min: sec 

Early 


HOLD 

Amount of time that airplane must delay 
before starting descent if crossing-time 
constraint at metering fix is to be 
satisfied, mln:sec 

Open-loop guidance in the form of desired altitude as a function of distance 
along the profile may also be computed by the pilot. This Is accomplished by 
keying a DME mileage Indication Into the display and pushing the "DME'VH" 
key. The desired altitude corresponding to that distance will then be com- 
puted and displayed to the pilot. This Information may then be used as 


guidance throughout the descent* 
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Figure 1. - Vertical plane geometry of computed descent path. 
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Figure 2. - Flight management descent algorithm logic flow. 
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Figure 6. - Modeled vertical speed at sea level for the DC-10 airplane. 
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Figure 7. - Flight data and generic parabolic model with c^ 
Mach number descents (idle thrust) for the DC 







220 knot 
airspeed 



Mach number 


Figure 8. - c. coefficient model for constant Mach number descents 
for the DC- 10 airplane. 
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Figure 9. - Modeled and actual true airspeed reduction in level flight 
(idle thrust) for the DC-10 airplane. 
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Figure 10. - Vertical speed at sea-level as a function of gross weight for 
the DC-10 airplane. 





Figure 11. - Progratnmabl e Descent Calculator. 
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